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The RHIC Beam Energy Scan focuses on the study of the QCD phase diagram — temperature (T) vs. baryon 
chemical potential (^s). The aim is to verify some predictions from QCD: that a cross-over occurs at fig = 0, 
that there exists a first-order phase transition at large fig and a critical point at an intermediate fig- The 
' spectra and ratios of produced particles can be used to extract T and fig in different energies and system sizes. 

' The STAR experiment has collected data for Au-|-Au collisions at ^sjvAf = 7.7 GeV, 11.5 GeV, and 39 GeV 

in the year 2010. We present midrapidity pT spectra, rapidity density, average transverse mass, and particle 
^1^, ratios for identified hadrons from the STAR experiment. The centrality and transverse momentum dependence 

O 1 of the particle yields and ratios are compared to existing data at lower and higher beam energies. The chemical 

■ and kinetic freeze-out conditions are extracted from the ratios and particle spectra. 

I I 1. Introduction 

>< 

The RHIC beam energy scan (BES) is devoted to exploring the QCD phase diagram [Il-i]. The QCD phase 
diagram represents the variation of temperature T vs. baryon chemical potential fj,B- At low temperatures, 
the relevant degrees of freedom are expected to be hadronic but at high temperatures, the quarks and gluons 
are the relevant degrees of freedom of the system. In the QCD phase diagram, at /Lte ^ 0, the transition from 
hadronic gas to quark gluon plasma (QGP) is expected to be a crossover At large /zb, it is expected to 
be a first order phase transition Q. The point where the first order phase transition line ends is called the 
QCD critical point which is of current interest for many heavy- ion experiments [H, The main objectives 

of RHIC BES are to search for the possible QCD phase boundary and to search for the possible QCD critical 
point in the QCD phase diagram. In addition, it will be interesting to see how various observables, suggesting 
QGP formation at top RHIC energies, behave as the collision energy is decreased. This will help to locate the 
beam energy where there is no QGP formation or the energy representing "turn off" of the QGP signatures. A 
^\ few potential observables for exploring the QCD phase diagram are: constituent-quark scaling of elliptic flow, 
[ parity-violation in strong interactions, wiggle shape of directed flow of protons, freeze-out eccentricity, K/tt 
ratio fluctuations, and higher moments of net-protons pj. 
" . I The QCD phase diagram can be accessed by varying temperature T and baryonic chemical potential /is- 

Experimentally this can be achieved by varying the colliding beam energy. A T-fi^ space point in the QCD 
phase diagram can be obtained from the spectra and ratios of the produced hadrons. Once the T-fi^ space 
point is obtained, one can study various signatures, mentioned above, for evidence of the possible QCD phase 

■ boundary and QCD critical point. In the year 2010, the STAR experiment collected data for Au-|-Au collisions 
at energies ^/SNN — 7.7, 11.5, and 39 GeV, as part of the proposed energies for the BES program, in addition 
to usual energies of 62.4 and 200 GeV. 



> 



2. Data Analysis 

The data presented here are for Au-|-Au collisions at y^SNN — 7.7, 11.5, and 39 GeV taken in 2010 by the 
STAR experiment Q. The total events analyzed for the results presented here are about 4 M, 8 M, and 10 M, 
respectively for ^/snn = 7.7, 11.5, and 39 GeV. The trigger selection is done by using the Beam Beam Counter 
(BBC) and the Vertex Position Detector (VPD) Q. The main subsystem used for particle identification is the 
Time Projection Chamber (TPC) The particle identification is enhanced to higher pT with the inclusion 
of the full barrel Time Of Flight (TOE) detector [l^]- The centrality selection is done using the uncorrected 
charged track multiplicity measured event-wise in the TPC within \ri\ < 0.5. The centrality classes represent 
the fractions of this multiplicity distribution. The average number of participating nucleons ((iVpart)) and 
collisions ((-/Vcoii)) are obtained by comparing the multiplicity distribution with that from Glauber Monte Carlo 
simulation The raw yields are extracted at low pT using ionization energy loss {dE/dx) from TPC, and 
using TOE information at higher pT- The raw spectra are corrected for the detector acceptance and tracking 
efficiency. These correction factors are obtained together by embedding the tracks simulated using the GEANT 
model of the STAR detector into real events at the raw data level [ill ■ The results are presented for midrapidity 
\y\ < 0.1 region. 
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Figure 1: Transverse momentum spectra for charged pious tt^ (left panel) and charged kaons K"^ (right panel) for various 
collision ceutralities at midrapidity {\y\ < 0.1) in Au+Au collisions at y^sjviv ~ 39 GeV. Errors shown are statistical. 



3. Results and Discussions 

3.1. Transverse Momentum Spectra 

Figure [T] shows invariant yields d^N/{2TTpTdydpT) vs. transverse momentum for tt^ (left panel) (right 
panel) and for various collision centralities at midrapidity (|y| < 0.1) in Au+Au collisions at ^/snn — 39 GeV. 
The curves represent fits to the spectra: pion spectra are fitted with Bose-Einstein function and the kaon spectra 
with TTiT-exponential. The dN/dy and (pt) or (mx) values are obtained using the data in the measured pr ranges 
and extrapolating using a fit function (Bose-Einstein for tt and mr-exponential for K) for the unmeasured pt 
ranges. The contribution to the yields from extrapolation to the total yield is about 20-30%. The pion spectra 
presented here are corrected for weak-decay feed-down and muon contamination using STAR HIJING+GEANT, 
as was done for previous STAR results [l2|. The total contribution from weak-decay feed-down and muon 
contamination is about 17% at low px and becomes almost negligible (~ 1%) around pr — 1-4 GeV/c. 

3.2. Centrality Dependence of Yields and Average Transverse Momentum 

Figure [5] shows the centrality dependence of dN/ dy normalized by number of participating nucleon pairs 
(A^part)/2 (top panels) and {pt) (bottom panels) for tt^ (left) and for (right) in Au+Au collisions at 
^ysNN = 7.7 and 39 GeV. New measurements from BES data are compared to previously published results at 
top RHIC energies The dN/dy per participating nucleon pair for pious is almost constant as a function 
of collision centrality at ^/snn = 7.7 GeV, while it increases with A'pait for higher energies and for kaons. The 
(pt) increases with collision centrality for both pions and kaons at all energies, indicating that the average 
collective velocity in the radial direction increases with collision centrality. 



3.3. Energy Dependence of Yields and Average Transverse Mass 

Figure [3] shows the energy dependence of dN/dy normalized by (iVpart)/2 (top panels) and (mr) — ni (bottom 
panels) as a function of ^/snn for tt^ (left) and for (right) . The star symbols represent the data from BES 
energies and are compared with previous measurements [12l - ll4| . The results are shown for the central collisions 
and for midrapidity regions. The STAR BES results are consistent with the published beam energy dependence 
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Figure 2: Left: Variation of dN/dy divided by participating nucleon pairs (top panel) and (px) (bottom panel) with 
collision centrality for tt*. Results from 7.7 and 39 GeV BES data are compared with those from previous STAR 
measurements [1^. Errors are statistical and systematic added in quadrature. Right: Similar measurements for . 
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Figure 3: Left: dN/dy divided by participating nucleon pair (top panel) and (mr) — m (bottom panel), plotted as a 
function of beam energy for vr^. Results from BES data are compared with those from previous measurements 
Errors are statistical and systematic added in quadrature. Right: Similar measurements for K^. 
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Figure 4: Left: Correlation of /K^ ratio with p/p ratio for central collisions at midrapidity, shown for new measure- 
ments from BES data and for the published results from RHIC energies _12|. The curve represents a the power law fit 
to the data. Errors are statistical and systematic added in quadrature. Right: Energy dependence of /C^/tt* ratio for 
central collisions at midrapidity. New results from BES data are compared to previously published results at RHIC p^ . 
SPS and AGS p^ . Errors are statistical and systematic added in quadrature. Results are also compared with 
various theoretical model predictions [20l-;24[ ] . 



trend. The yields per participating nucleon pair increase with beam energy. The quantity (my) — m, where 
Tot = v^'t + tt? and m is the hadron mass, increases with beam energy for lower energies, becomes almost 
constant for the region covered by the BES data and then tends to increase towards the top RHIC energies. 
This is an interesting observation in a scenario where the system is in a thermodynamic state. In that case, 
(tot) — rn can be related to temperature of the system and dN/dy (cx \og{^SNN)) may represent entropy. 
Then this observation could reflect the signature of a first order phase transition as proposed in the Ref. [151]. 
However, other interpretations of the observed (tot) — to are possible [l6| . 



3.4. Particle Ratios 

The particle ratios provide useful information about the collision dynamics as these are fixed after the chemical 
freeze-out (see Sec. 3.5). Figure |4] (left panel) shows the correlation of K~ jK^ ratio with p/p ratio for new 
measurements from BES energies along with our previously published results ■ This could give information 
on how the kaon production is related to the net-baryon density which is reflected by the p/p ratio. The 
STAR -p and p yields are not feed-down corrected. At lower energies, the kaon production is dominated by the 
associated production which results in more production compared to . Also the p/p ratio is much less 
than unity at midrapidity, indicating that there is large baryon stopping at midrapidity at the lower energies. 
As we go towards higher energies, the pair production mechanism starts to dominate and the ratios tend to 
become closer to unity. The correlation between K~ jK^ (representing net-strange chemical potential, [is) and 
p/p (representing net-baryon chemical potential, /is) ratios seem to follow a power law behavior with a=0.2 
(represented by the curve). In a hadron gas, the relationship between /zg and /is depends on the temperature. 
In a particular case of T=190 MeV and \ib < 500 MeV, these potentials follow the relation /is = (1/3)/xb [l7l |. 

Recent theoretical calculations [l^ suggest that the maximum net-baryon density at freeze-out is attained 
at the lowest BES energy of ^/snn ~ 7.7 GeV. The maximum net-baryon density could also be related to the 
peak observed in the energy dependence of /tt^ ratio at around ^^snn ~ 7-8 GeV, as was observed by the 
NA49 experiment [l^l- This is sometimes referred to as the "horn". The K/n ratio is also important as it could 
suggest the strangeness enhancement in heavy- ion collisions with respect to the elementary collisions. Figure S] 
(right panel) shows the energy dependence of if^/7r* ratio for central collisions at midr apidity . The BES results 
are shown with star symbols and are compared to results from previous measurements [12h14| | . The BES results 
are in good agreement with the trend of energy dependence established by the published measurements. 

Figure |4] (right panel) also shows the predictions of energy dependence of K/n ratio from various theoretical 
model calculations. The energy dependence of /tt^ ratio has been interpreted using the Statistical Model 
of Early Stage (SMES) [l^ . The model predicts first order phase transition and the existence of mixed phase 
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Figure 5: Energy dependence of baryonic chemical potential (top-left), temperature (Teh, rkin) (top-right), and average 
flow velocity (bottom-left). Bottom-right: Chemical freeze-out temperature plotted as a function of baryonic chemical 
potential. New results from BES data are shown with star symbols while other results are from the Ref. and 
references therein. Curves represent the parametrization from Ref. [2?! ]. 



around beam energy of 7-8 GeV. The SHM or Statistical Hadronization Model [2U\ assumes that the strong 
interactions saturate the particle production matrix elements. This means that the yield of particles is controlled 
predominantly by the magnitude of the accessible phase space. The system is in chemical non-equilibrium for 
y/SNN < 7.6 GeV, while for higher energies, the over-saturation of chemical occupancies is observed. The 



Statistical Model [21[ assumes that the ratio of entropy to as a function of collision energy increases for 
mesons and decreases for baryons. Thus, a rapid change is expected at the crossing of the two curves, as 
the hadronic gas undergoes a transition from a baryon-dominated to a meson-dominated gas. The transition 
point is characterized by T=UO MeV, 410 MeV, and y/sWN=8.2 GeV. In the Thermal Model [H, the 

energy dependence of /tt^ is studied by including cr-meson, which is neglected in most of the models, and 
many higher mass resonances (m > 2 GeV/c^) into the resonance spectrum employed in the statistical model 
calculations. The hadronic non-equilibrium Kinetic model [2^ assumes that the surplus of strange particles is 
produced in secondary reactions of hadrons generated in nuclear collisions. Then the two important aspects are 
the available energy density and the lifetime of the fireball. It is suggested that these two aspects combine in 
such a way to show a sharp peak for the strangeness-to-entropy or K/tt ratio as a function of beam energy. In the 
Hadron Resonance Gas and Hagedorn model (HRG-HHagedorn) [2J], all hadrons as given in PDG with masses 
up to 2 GeVc^ are included. The unknown hadron resonances in this model are included through Hagedorn's 
formula for the density of states. The model assumes that the strangeness in the baryon sector decays to strange 
baryons and does not contribute to the kaon production. The energy dependence of K'^/n'^ ratio seems to be 
best explained using HRG-|-Hagedorn model. 



6 



Proceedings of the DPF-2011 Conference, Providence, RI, August 8-13, 2011 



3.5. Freeze-out Conditions 

The integrated yields of hadrons change only through the inelastic collisions. The point in time when these 
inelastic collisions stop is called the chemical freeze-out. After this stage, the particle ratios are frozen and the 
system can be described by the thermal equilibrium model [T2| . [25 !| which involves two main parameters Teh 
and /is. The Tch and /is are obtained by fitting the particle ratios with the thermal model. The point in time 
when the the elastic collisions among the particles cease is called the kinetic freeze-out. The system can be 
described by the blast- wave (BW) formulation involving two main parameters Tkin and average flow velocity 
((/?)) [T^. [26j. The Tkin and are obtained by simultaneously fitting the invariant yields for tt, K, and p with 
the blast- wave model. 

Figure [5] shows the variation of extracted freeze-out parameters for central collisions. Top-left plot shows the 
energy dependence of baryonic chemical potential. It shows that ^lb decreases as the beam energy increases. 
This is expected since there are fewer net-baryons at midrapidity at higher energies because of less stopping 
of baryons at midrapidity at higher energies. Top-right plot shows the kinetic freeze-out (open symbols) and 
chemical freeze-out (solid symbols) temperature plotted vs. ^/snn- The chemical freeze-out temperature 
increases with energy and saturates at the higher energies, and kinetic freeze-out temperature decreases with 
beam energy after ^/snn ~ 7.7 GeV. Bottom-left plot shows the average flow velocity plotted vs. y^s^vw- The 
increases with beam energy. Bottom- right plot shows the current picture of the phase diagram (Tch vs. /is). 
The new measurements from BES program now extend the /is range covered by the RHIC from 20-400 MeV. 

4. Summary 

In summary, bulk properties from the RHIC BES program are presented in this paper. The BES data are 
nicely fitted in the beam energy dependence trend established previously. The yields of identified hadrons 
increase with beam energy. The {tut) — m is almost constant for the BES energies. The (pt) increases with 
centrality reflecting collectivity increases with centrality. The net-baryon density plays an important role at the 
lower energies as seen in the energy dependence of K/tt ratio, and the correlation between K~ / and p/p 
ratios. The new measurements from the BES program extend the hb range covered by the RHIC from 20-400 
MeV. 
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